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Abstract 

In many systems in condensed matter physics and quantum field theory, first order 
phase transitions are initiated by the nucleation of bubbles of the stable phase. In homoge- 
neous nucleation theory the nucleation rate P can be written in the form of the Arrhenius 
law: P = Ae~'^''. Here TLc is the energy of the critical bubble, and the prefactor A can 
be expressed in terms of the determinant of the operator of fluctuations near the critical 
bubble state. In general it is not possible to find explicit expressions for A and Tic- If the 
difference r] between the energies of the stable and metastable vacua is small, the constant 
A can be determined within the leading approximation in r], which is an extension of the 
"thin wall approximation" . We have done this calculation for the case of a model with a 
real-valued order parameter in two dimensions. 

PACS numbers: 05.70.Fh, ll.lO.Kk, 64.60. Qb 
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1 Introduction 

The problem of the decay of the metastable false vacuum at first order phase transitions has 
attracted considerable interest due to its numerous relations with condensed matter physics 
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quantum fields 0, cosmology [Q, and black hole theory 0]. In hanger's theory of homogeneous 
nucleation the false vacuum decay is associated with the spontaneous nucleation of a 

critical bubble of a stable phase in a metastable surrounding. In the context of quantum field 
theory, the nucleation theory was developed by Voloshin et al. 0, and Callan and Coleman 
0. The quantity of main interest is the nucleation rate / per time and volume. 
The nucleation rate 

ZTT 

is a product of the static part F and the so-called kinetic prefactor k, which depends on 
the detailed non-equihbrium dynamics of the model, see @, Most important is the static 
nucleation rate F, which is equal to twice the imaginary part of the free energy density of the 
metastable phase. In this article we study the static part F. 

In the homogeneous nucleation theory it has the form of the Arrhenius law: 

F = ^exp(-K), (2) 

where Tic is the energy of the critical bubble. The prefactor A is determined by fluctuations 
near the critical bubble state and can be expressed in terms of the functional determinant of 
the fluctuation operator 0, |^ . 

In the general case, it is not possible to find the explicit critical bubble solution of the field 
equations analytically. However, the problem becomes asymptotically solvable, if the decaying 
metastable state is close enough in energy to the stable one, i.e. if the energy density difference 
rj between the metastable and stable vacua is small. The leading approximation in this small 
parameter is usually called the "thin wall approximation" |TD|, since at — the critical 



bubble radius goes to infinity and becomes much larger than the thickness of the bubble wall. 

In the thin wall approximation, the critical bubble energy Tic can be easily obtained from 
hanger's nucleation theory. It turns out to be much more difficult to find explicitly the prefactor 
A in Eq. (j^). This problem, which is important for applications of nucleation theory, has been 
extensively studied in different models. 



A remarkable result on this subject was obtained by Voloshin W^. He considered scalar 
field theory in 2 dimensions with a potential f/(0) of the type shown in Fig. 1. 

Voloshin claimed that in the limit — > the nucleation rate F in such a model can be 
described by the simple universal formula 

^''] ■ (3) 

Here a is the surface tension of the wall between the stable and metastable vacua in the limit 




7] ^ Q. Thus, according to [|TT|, in this limit the nucleation rate F is determined by two well 



defined macroscopic parameters rj and cr. Another claim of |]TT| is that there are no corrections 



to Eq. proportional to powers of the dimensionless parameter rj/a'^. Voloshin arrived at 
these conclusions by an analysis performed in the thin wall approximation. He replaced the 
original scalar field theory by an effective geometrical one, which describes only fluctuations 
of the critical bubble shape. This approach implies that all other fluctuations of the original 




Figure 1: The potential U with the false 



and true 



vacuum 



scalar field could be properly accounted for by the correct choice of the macroscopic parameters 
r] and a. 



Recently an analytical method was developed [12|, which allows one to study nucleation in 



the scalar field model beyond the thin wall approximation. In [12] this method was used to calcu- 
late the nucleation rate for the first order phase transition in the three-dimensional Ginzburg- 
Landau model. In the present paper we apply the same approach to the two-dimensional 
case. We calculate the nucleation rate beyond the thin wall approximation and verify directly 
Voloshin's claim (^. 

Nucleation theory in two-dimensional scalar field theory has also been studied by Kiselev 
and Selivanov Strumia and Tetradis |T^, and other authors. In these articles, however, 



different renormalization schemes have been used and F has not been expressed in terms of 
macroscopic parameters t] and a. This makes it difficult to compare their results with the ones 
discussed in this article. 



In the articles [T^, [T^ the nucleation rate was calculated in the two-dimensional Ising model 
in a small magnetic field for arbitrary anisotropics. If Voloshin's result @ is universal, it should 



be applicable as well to the Ising model in the critical region. Indeed, the results of [16 



rewritten in terms of f] and a, are in a very good agreement with Eq. (|^). The exponent factors 
are the same, and the prefactors differ only by the number 7r^/9 ~ 1.0966, which is very close 
to unity. This small discrepancy increased our interest in the subject of the present study. 

2 Model and notations 



We consider the two-dimensional asymmetric Ginzburg-Landau model defined by the Hamilto- 
nian: 



d X 



{d,<P{x)f + U{<P{x)) 



(4) 



where is the continuous one-component order parameter, and the potential f/(</)) depicted 
in Fig. 1 is given by 

Ui<p)=U,i<P) + ^{<P-v) + Uo. (5) 
Here Us{(j)) denotes the symmetric part of the potential: 



(6) 



The potential U (0) has a metastable minimum (false vacuum) at (p = (p^ and a stable one (true 
vacuum) at = The constant term Uq in Eq. (|^) is chosen to ensure t/(0+) = 0. 
The partition function is given by the functional integral 



D(j) exp [-ni(j))] . 



The temperature has been absorbed into Ti. 

It is convenient to define the mass m and the inverse coupling parameter /3 by 



92 



m 



■Usicp) 



3m 



and to introduce dimensionless quantities 



9 

— r/0, ^{x) = 



0± ^ ^ 



2 m'* ■ ■ V V 

In dimensionless variables the Hamiltonian and partition function take the form 



d X 



where 



and 



D(p{x) exp —/3H{(p) 



(7) 



(8) 

(9) 

(10) 

(11) 
(12) 



3 The critical bubble solution 

The uniform solutions of the field equation 

sn/Sip{x) = (13) 

are the stable if- and false (metastable) vacua given by 

¥'± = ±l-|T|-^ + 0(r^^). (14) 

The critical bubble ipb{x) is the non-uniform radially symmetric solution of Eq. ([T3| ) approaching 
the false vacuum at infinity. That is, 

- ^ - + 2Mvl -l) + U = 0, lim Mr) = (15) 




Figure 2: Profile of tlie critical bubble 

where f = ^5^7^- The profile of the critical bubble solution is shown schematically in Fig. 2. 
If fj is small, the thin wall centered at f = i? divides regions of false and stable vacua outside 
and inside the bubble, respectively. 

Equation ([15|) can not be solved explicitly. Following the approach introduced in |]l2l we 



shall construct the solution by expansion in powers of 7]. Introducing the new independent 
variable C,'- 

^ = f-R, (16) 

we expand R and ipbiO 

R = ^+00 + 01^7 + 02^^ + 0(7]^), (17) 

MO = MO + MOv + MOv' + oiv'). (i8) 

After substitution of Eqs. (p!6|-[18D into (pTSj) one obtains perturbatively in fj: 

12 , , 

a_i = 2' ^0 ^ ^' ^ ~9' = 0) (19) 

MO = tanh^, ¥?i(0 = ~, 

'^2(0 = -— ^-27|lOe-16ecosh(20-2ecosh(40 + 
24 cosh 4 I 

21n 12cosh«] m + 8-h(2«) + slum)] - 2, jit nanh^}. 



The bubble energy E = H [iphix)] can be written as 

£ = ./de(fl + o(^)'. (20) 

-R 



Substitution of Eq. ([ID into (^) yields 



E 



2n 



1 



-+V , ,r 



19 TT^ 



+ Oif] 



(21) 



It is the basic principle of homogeneous nucleation theory that the decay of the metastable 
vacuum occurs through nucleation of the critical bubble. Callan and Coleman expressed the 
nucleation rate F of the metastable vacuum in terms of functional determinants §]. In our 
notation their result takes the form 



PE 



2- 



= exp(-/3E + 5) 



(22) 



Here T = AV/rn? is the dimensionless nucleation rate, and the entropy S associated with the 
critical bubble is given by 

det'M 



exp S 



-1/2 



(23) 



det M(o) 

where M and M^^^ are the fluctuation operators near the bubble v^fe(5;) and the metastable 
uniform vacuum respectively: 



M 



-d' + 6 [^,{r)Y - 2 
-9^ + 6(^^-2. 



(24) 
(25) 



The operator M has two zero modes proportional to ^h{x)i fJ' = ^,"2, and one negative mode 
with the eigenvalue 

Ao = -4f. (26) 

The notation det' implies that the three above mentioned modes are omitted in the correspond- 
ing determinant. After substitution of Eqs. (|2ll) and (^), equation (^) simplifies to 



P 



— — exp 

Qr]^ \ 3r] 



^""^ + 3] il + 0{f,)). 



(27) 



In the subsequent sections we shall calculate the small rj expansion for the critical bubble 
entropy ( pSf ) to the order 0{fj^). 



4 The bubble entropy 

The spectrum of the fluctuation operator M can be determined in the form of a perturbative 



expansion in powers of the parameter rj as in |12|. This is achieved in the following way. 



Introducing the angular momentum quantum number G Z in two dimensions, the radial 
Schrodinger operators corresponding to M and M^^^ are 



- -|5-^;^ + ^H-6MOf-2, (28) 

(f Id /i^ 

df'^ f df 



= -^-54 + ^ + 6/;-2. (29) 



Shifting the coordinate from f to and making use of the Laurent series ([T7|) for R, the 
eigenvalue problem for if^ can be treated perturbatively in fj. The lowest order leads to the 
exactly solvable Poschl- Teller operator 

-^-Qsech^^ + A, (30) 

which has discrete eigenvalues and 3 and a continuum above 4. In second order one finds a 
band around 0: 

Xo, = 4f{f,'-1) + O{f), (31) 



a band around 3: 



and a continuum 



A3^ = 3 + 4?7"(/i' + const.) + 0(r7^), (32) 



Xk^ = e + 6^1-2 + 477 V' + 0{f]^) (33) 

with e R. This spectrum includes the negative mode Aqo and the two zero modes Ao,±i. 

The sum over /i and the integration over k produce ultraviolet divergencies in S. We treat 
these by means of dimensional regularization in = 2 — e dimensions. As there appear volume 
integrals in intermediate steps of the calculation, the extra dimensions are equipped with a 
finite extent L and periodic boundary conditions. The parameter L must cancel out in finite 
results. 

The finite part of the regularized entropy can be conveniently evaluated with the help of 
zeta-function techniques Jl^ . The operator-zeta function appropriate for our case is defined by 



1 



00 



Cuiz) = / (Tr 'e'''' - Tr e-**^'"' ) (34) 

1 [z) Jo ^ ' 

for Re^; > 1 and analytical continuation to other values of z. The integrand contains the 
heat kernels exp(— tM) and exp(— tM*^°^). For positive t there is an asymptotic expansion, the 
so-called Seeley expansion, which is of the form 

00 

Tr (e-**^ - e-**^*"^) = (47rt)-'^/2 g _ (35) 

n=l 



Following one obtains 

1 A n. r9 1 

0{e). (36) 



^=^;fcM(0) + ^ 
2 az Stt 



^ + ln47r + r'(l^ 



This expression displays the divergence as a simple pole in e. The derivative of the zeta-function 
is a finite quantity. The first Seeley coefficient is given by 



= _6 {[Mr)? - vl) =L-^^ + 0{v), (37) 

J 7] 

where 



The zeta-function is decomposed into a contribution from the band near zero and the rest, 



where 



Co(;^) 



Cm(^) = Co(^) + Ci(^), 



1 



V{z) 



poo 



Correspondingly, the entropy is decomposed as 



S = So + Si + L-'— 
477 



- + ln47r + r(l) 
e 



+ 0{fi) + 0{e). 



Consider the part 

For a general spectrum of the type 



= a{ii + b){ii + c) 



(39) 
(40) 

(41) 
(42) 
(43) 



the logarithm of the zeta-function regularized determinant is given by 

- -^Co(O) = -21nr(6 + l) -21nr(c+l) - (6 + c) In a + ln(6c) + 2 ln(27r), (44) 
az 

which can be derived with the help of Riemann's and Hurwitz's zeta-functions. Setting a = 
Arf, b = 0,c = 2 one finds 

This piece will therefore contribute the factor 

877^ 



TT 



(45) 



TT 



{1+0(77^)} 



(46) 



to the prefactor A of F. 

The remaining part 5*1 of the entropy is calculated with the help of methods from quantum 
mechanical scattering theory. The heat kernels can be represented as 

Kt{M) = Tr (e-*^ - e"*^^"') = - / ^ e"^* Tr [(A - M)-' - (A - M^)"!] , (47) 

c 

where the integration path C in the complex plane is shown in Fig. 3 
Decomposed into the angular momentum sums this reads 



with 



A{\ix)^T.[{X-H,r'-{\-Hf)-']. 



(48) 



(49) 



We obtained an exact representation for A{\,ii). To describe it some notations are necessary. 



c 

^ ► 



Figure 3: Integration path C in the complex A-plane 



Let fiif, A), gi{r, X), i = 1, 2 be the solutions of the linear ordinary differential equation 

= A ^(f) (50) 

determined by their asymptotics: 



/i(r, A) Kf,{qr), and /2(r, A) If^{qr), at r ^ oo, 

gi{f,X) f), and 5f2(r, A) ^ J/,(g_ f), at f ^ 0. 



(51) 
(52) 



Here K^{z) and I^{z) are modified Bessel functions, and the parameters q and g_ are defined 
as 



A-2 



1/2 



A-2 



1/2 



(53) 



Since the second order equation (^) has two linearly independent solutions, there is a linear 
dependence between the functions (^^(f. A) and /j(f. A): 



The function A{X, fi) can be expressed explicitly in terms of the coefficient a22(A): 



(54) 



(A + 2-6<y92)(A + 2-6(^^ 



+ 



d\n 0:22 (A) 

dx ■ 



(55) 



This representation for the trace of resolvent operators is exact. However, equation 
can not be solved in closed form for arbitrary fj. So we have to consider the small-?] expansion 
for ^22 (A). We have obtained two terms of this expansion by use of a perturbation theoretical 
analysis of the scattering problem (^Oj-Q). Omitting the details, the logarithmic derivative of 
the matrix element 022 (A) up to quadratic terms in fj takes the form: 



(ilna22(A) 
dX 



(4 + p2 _ xy/2 



2p2 



(A -4)' A-4 A-3-J92 A-p2 



+ Oifi'). (56) 



Here p is the angular momentum parameter defined as p = 2r]fi ^ fi/ R. 



Substitution of Eqs. (H) and (|T4D into (H) yields: 



A{X,fi) 



1 



2p2 



(4 + p2_^)i/2 [(A -4)^ A -4 

2 \P\ , ^..2 



+ 



A-3 



+ 



A — 



(A-4r 



0(r/^). 



(57) 



This function has simple poles at A = p^, A = 3 + and a square root branching at A = 4 + 
It is analytic at A = 4. 

With the help of this expression we could evaluate the integral representation for Kt{M) 
and the related zeta-function. The details of this lengthy analysis will not be presented here. 
The /i-summations have been done by means of Poisson's summation formula. Separating 
the contribution of the band near zero, which has been treated above, the final result for the 
remaining entropy is 

(58) 



This is the central result of this section. 



5 Decay rate 

In d = 2 — e dimensions the bubble energy is associated with an additional factor and the 
regularized dimensionless decay rate F is given by 



_ ■ exp z — 

6 7]^ \ 3r] 

In terms of the dimensionless quartic coupling 



+ s]{i+om 



u 



9 



m 



A-d 



the parameter (3 is equal to 



R 3 _ 

p = —m 
u 



Using the results (0,|^), we can write the nucleation rate as 

"2 



(2L-n^exp 
nu 



(2L)- + L-^-^ 



+ ln47r + r'(ll 



+ Si + 0{fi) + 0{e) 



(59) 

(60) 
(61) 

(62) 



The entropy contains an UV-divergent term, represented by a pole in e. After renormalization 
of the parameters of the model according to the usual prescriptions, the divergencies as well as 
the spurious L-dependence should disappear in the hmit e ^ 0. 

For the renormalization of the model parameters we use the same scheme as in [|l^, [T^. A 
straightforward calculation yields the relation between the bare and renormalized dimensionless 
couplings and masses on the one-loop level: 



u 



) 1 



- + ln47r + r(l) + ^ 

e 4 



0{e) 



(63) 



m = <i 1 + 



Mr 
47r 



- + ln47r + r(l) + ^ 

e 4 



+ 0{e) + 0{ul) 



(64) 



The asymmetry parameter r^o is renormalized as follows. On tree level the difference between 
the minima of the potential 

At/ = - = r/o + 0{r,l) (65) 

is equal to ?7o for small asymmetries. Let Ues be the full effective potential given by 



(fx f/cfr(0) 



for 



const., 



(66) 



where r[(/)] is the generating functional of one-particle irreducible vertex functions. Then we 
define the renormalized asymmetry parameter rj through 



r] = Af/eff = t/efr((0) + ) - f/cff((0)-). 



It is related to r/o by 



r/ = ^(0) + O(r^o'), 



(67) 



(68) 



where (0) is the expectation value of the field at rjo = 0. From a one-loop calculation we get 



Vo 



- +ln47r + r(l) 



+ 0(s) + 0(4)}. 



(69) 



Expressing the unrenormalized parameters in terms of their renormalized counterparts, the 
divergencies cancel indeed and in the limit e = we obtain 



with 



F ^ -.,+4.^{l + ^g-51n2)+0(4)}+0(.) 



27r 



rrif. 



(70) 

(71) 
(72) 



To make contact with Voloshin's proposal we have to express this result in terms of the 
interface tension a. To this end we calculated a along the lines of [17, 18 1, but now for d = 2. 
Leaving out the details here, we obtained 



a 



Ur 



This implies 

Our final result is therefore 

r 



TTO" 



{l + 0{ul)} + 0{v). 



2tt 



exp 



[l + 0{ul)} + 0{^) 



(73) 



(74) 



(75) 



This is in perfect agreement with Voloshin's result. 

Studies of the nucleation rate in the two dimensional Ising model have been made with the 
Monte Carlo method, see e.g. |20|, We would like to add a remark on this case. Let 



Z = J2exp\KY^ S,Sj + Hj2sA (76) 

{S} [ <ij> i J 

be the partition function for the two-dimensional Ising model on a square lattice with lattice 
spacing 1. The critical coupling is 

K,= Uyi{V2 + 1). (77) 

The quantities appearing in Eq. ([75| ) can be related exactly to K and H in the critical region. 
The interface tension is given by pT|] 



(T = 2ir + lntanhir (78) 

and the asymmetry parameter is 

7] = 2MH, (79) 

where 

M = (l - [sinh 2K] (80) 



is the zero field magnetization P^ |. 

The kinetic prefactor k, mentioned in the introduction, cannot be calculated with static 
methods, because it depends on the non-equilibrium dynamics. For dynamics that can be 
described by a Fokker-Planck equation, it is, however, expected to be proportional to the 
negative eigenvalue |Ao| and contributes two additional powers of the magnetic field |p, 

6 Summary 

Our semiclassical calculation of the nucleation rate F in the two-dimensional Landau- Ginzburg 
0^-model confirms Voloshin's result (^, which was derived in the thin wall approximation. In 
particular, we confirm the prefactor value A = ri/{2'K) first obtained by Kiselev and Selivanov 
1l3| , and Voloshin [0 . 



This value differs from that obtained for the two dimensional critical Ising model |T^, |T| 
by the numerical factor vr^/Q ~ 1.0966. We suppose that this small discrepancy is the result of 
approximations used in |jl5|, |l6l, and the prefactor value A = ?7/(27r) is universal. 
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